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ABSTRACT 

There  were  three  directions  to  this  study-  theoretical  modeling,  collection  and  systematization  of 
experimental  data,  and  analysis  of  long-term  variations  in  the  Arctic  Basin.  The  work  with  the 
experimental  data  included: 

1.  The  creation  of  a  new  version  of  the  Arctic  Basin  Data  Base  by  combining  the  WOA-94, 
MOODS,  and  the  historical  Russian  data.  The  data  base  will  be  used  for  studying  the 
evolution  of  large-scale  thermohaline  anomalies,  as  well  as  further  development  of  the 
theoretical  model. 

2.  The  collection  and  preliminary  analysis  of  the  archive  data  on  the  water  exchange  through 
Fram  Strait  and  the  runoff  of  the  Siberian  rivers  as  important  factors  influencing  the 
oceanographic  and  ice  conditions  in  the  Arctic  Ocean.  (These  data  will  be  used  later  for 
modeling  of  the  Arctic  Ocean  response  to  external  impacts). 

3.  A  comparison  of  the  heat  content  variability  in  the  Atlantic  layer  in  different 
oceanographic  regions  outlined  in  the  report  "Climatic  variability  of  the  Arctic  ocean" 
(Final  report,  MSIC,  1996). 

4.  A  comparison  of  salinity  variations  in  the  upper  ocean  layers  in  these  regions. 

The  aim  of  the  theoretical  work  was  to  improve  the  numerical  scheme  for  modeling  the 
Atlantic  water  inflow  into  the  Arctic  Basin,  and  was  performed  by  N.G.  Iakovlev.  The  results  of 
the  modeling  are  described  in  Section  I  of  the  Report. 

A  description  of  the  new  version  of  the  Arctic  Basin  Data  Base  is  presented  in  Section  II. 
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I.  ON  NUMERICAL  MODELING  OF  THE  ARCTIC  OCEAN  RESPONSE 

TO  EXTERNAL  FORCING 


1. 1.  The  development  of  the  new  numerical  scheme 

1.1.1.  Rationale 

One  of  the  prominent  features  of  the  water  circulation  in  the  Arctic  Ocean  is  the  existence  of 
narrow  coastal  jets,  one  of  which  is  the  West  Spitsbergen  current  supplying  the  Arctic  with  warm 
Atlantic  water.  In  the  previous  version  of  our  general  ocean  circulation  model,  this  current  was 
almost  absent  because  of  comparatively  low  spatial  resolution  and  nonslip  boundary  conditions 
for  the  velocity  components  at  the  side  boundary,  and  at  the  bottom.  The  width  of  the  Fram  Strait 
is  about  400  km,  i.e.  ~  four  grid  steps.  To  allow  for  the  Atlantic  water  inflow  we  had  to  specily  a 
model  of  the  bottom  bathymetry  deeper  than  that  in  reality,  thus  making  it  necessary  for  us  to 
develop  a  new  numerical  scheme. 

The  first,  direct,  way  to  improve  the  simulation  of  the  coastal  jets  is  to  refine  both  the 
horizontal  and  vertical  spatial  discretization.  This  demands  no  changes  in  computer  programming 
and,  at  first  glance,  may  seem  very  attractive.  However,  even  when  we  choose  a  grid  size  three 
times  smaller,  it  barely  allows  us  to  get  results  on  the  circulation  sensitivity  in  reasonable  “real” 
time.  Moreover,  the  nonslip  boundary  condition  for  the  velocity  is  acceptable  only  when  the  grid 
size  is  much  less  than  the  internal  Rossby  radii  -  for  the  Arctic,  it  is  about  10  km,  i.e.,  10  times 
less  than  the  grid  size  used  in  the  model. 

The  second  way  is  to  design  a  special  coordinate  system  refining  the  areas  of  interest.  The 
simplest  versions  of  this  approach  have  been  proposed  by  S.  Hakkinen  (USA)  and  by  G. 
Semjonov  (AARI,  Russia).  The  approach  is  fine,  but,  unfortunately,  the  area  considered  changes 
in  time. 

We  decided  to  change  the  boundary  conditions  and  assume  conditions  with  free  slip  at  the 
side  boundary  and  quadratic  friction  at  the  bottom.  This  approach  does  not  contradict  the  first 
and  second  approaches,  but  may  even  complement  them,  if  necessary. 


1.1.2.  New  formulation  of  the  general  circulation  model 

The  basic  features  of  the  new  model  are  the  same  as  in  the  previous  version,  but  the  boundary 
conditions  for  the  horizontal  components  of  the  velocity  are  different. 

Let  us  consider  the  system  of  equations  of  large-scale  ocean  dynamics  under  commonly  used 
assumptions  in  a  3D  area  Q  bounded  by  plane  z  =  0  and  sufficiently  smooth  bottom  topography 
z=H(x,  y).  To  simplify  the  presentation,  we  use  the  Cartesian  coordinate  system  (x,  y,  z). 

du  IdP 

-—  +  Lm  -  /v  = - —  -t-  Du 

ot  Pq  ox 

— — Lv  +  lu  = - -I-  Dv 

dt  po  dy 

dP 
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(1) 


du  3v  3w  - 


dt 


+  LT=DtT 


^+LS  =  DsS 
dt  ^ 

P=P(r,s) 

Here  u,  v,  w  are  the  components  of  the  velocity  u  along  the  axes  x,  y,  z  respectively,  P  is  the 
pressure,  T,  S  are  the  temperature  and  the  salinity,  p  is  the  density  anomaly  compared  with  the 
constant  ,  /  is  the  Coriolis  parameter,  g  is  the  gravity  acceleration,  and  t  is  time.  L,  D,  Dt,  Ds 
are  the  operators  of  advection,  turbulent  mixing  of  momentum,  and  turbulent  diffusion  of  heat 
and  salt,  respectively. 

We  specify  the  advection  as 

ax  ay  az 

and  turbulent  mixing  and  diffusion  as 

9  3(|)  B  5(1)  3  9(1) 

^  dz  dz  dx^dx  dy^dy 

9  9(1)  9  9(1)  9  9(1) 

Let  the  boundary  consist  of  two  parts:  F  =  Fq  UF,,  where  F^  is  a  “liquid”  boundary  (passages 
and  rivers).  As  for  the  first  part  of  the  boundary,  let  Fq  =  F^Uf^o*  where  Qq  is  an  intersection 
of  Cl  with  the  plane  z  =  0,  and  Fg  is  a  “solid”  boundary  representing  the  land. 

Let  the  boundary  conditions  for  momentum  be  as  follows: 


(x) 


( y)  T- 

vpo^  =  x^^^  onFp, 


u„=0 


on  F, 


S  ’ 


(2) 


o 

w  =  —  on  Qq  , 
dt 

u  =  Uh  ,  V  =  V/,  on  Fj  . 

For  temperature  and  salinity  we  specify  the  following  boundary  conditions: 

^  dS  ^  ^ 

T=Th  ,  5  =  on  Fi  . 

Here  ^(x,  y,t)  =  — ^  Pix,  y,0,t)  is  the  sea  surface  elevation,  n  is  the  outer  normal  vector,  u„  is 
Po^ 

the  velocity  component  along  the  vector  n  ,  -:^  =  |i.-^cos(/i,x)  +  |i— cos(n,y)  +  v— cos(n,z) 

dN  dx  dy  dz 

is  the  conormal  derivative  corresponding  to  the  momentum  turbulent  mixing.  Derivatives 


dNi 
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and  - -  are  specified  in  the  same  way.  Functions  ,  v^,  7^,  5^,  Q7,  and 

BNs 

are  assumed  to  be  known. 

The  initial  conditions  at  t  =  0  are  as  follows: 

u  =  ,v  =  v^^\  T=T^^K  5  =  5^^^  and  C  = 

dll 

Note  1.  In  general,  the  linearized  kinematic  condition  w  =  on  z=0  is  valid  for  (^«H . 

at 

This  condition  is  not  satisfied  in  the  vicinity  of  F^finQ.  However,  in  a  discrete  model,  after 
applying  the  finite  element  technique,  we  find  min//(x,y)=  10m  ,  while  the  scale  of  level 
variations  is  about  Im. 


To  develop  a  finite  element  scheme  we  pose  a  projection  form  of  the  problem  (l)-(2),  where 

functions  K],K2,TtT^'^s  are  equal  to  zero  on  Fj  ,  and  functions  itj,  K4,  Tip  eWjiCl): 

du  I  dP  r  T  . 

(—,7C,)-B(u,U,Jt,)-(/v,7C,)  =  -  —  7C,)-[M,7t,J+/,  , 

9v  1  BP  r  ■] 

(— ,7t2)  -  fi(v,U,7t2)  +  (/m,K2)  =  “■^^■^’"2)  -  /2  - 

Bp 

(— ,7Cj)  =  (pg,;i5)  , 

Bz 

,  BK4  .  ,  BjIa  .  ,  Bltd .  . 

-  («,-—)  -  (V,-^)  -  (w,-^)  =  fs  , 

Bx  By  Bz 

BT 

(—,%)- 5(7, 11,117.  )  =  -[r,  Tty- +/7-  , 

BS  r  I 

(,—,Tts)-B{S,\i,its)  =  -[S,t^s\s+fs  ’ 

(p,7tp)  =  (p(7,S),jrp)  . 

The  following  denotations  are  used  here: 

(a,5)  =  J  abd^  , 
n 

r(  Bb  Bb  Bb^ 

B(a,u,b)  =  \\  au- — hav— +  flw—  dQ  , 

qV 

r  ,1  rf  BaBb  BaBb  BaBb^,^ 

and  [a,b]^  is  specified  in  similar  way.  The  right-hand  sides  fi,  fz,  fs,  fr,  and/s  are  generated  by 
boundary  conditions. 

AS  in  the  previous  version  of  the  model,  we  choose  right  triangular  prisms  to  be  the  finite 
elements.  Let  the  bases  of  the  prisms  lie  on  the  planes  z  =  z^,  k  =  1,...,Kq  ,  and  the  vertexes  of 
the  prisms  belong  to  the  set  of  points  Piix,y),i  -  .  Let  each  index  i  correspond  to  the 
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integer  number  KL(i)  <  Kq  ,  hence  will  be  treated  as  a  depth  at  point  P{i) .  On  the  finite 

elements,  we  specify  the  basis  functions  71,;^  =(Pi(x,y)\^ i^iz)  ,  where  (p/CJt.y)  are  2D  finite 
linear  piecewise  functions  defined  on  triangular  finite  elements  and  tl/jt(z)  are  ID  finite  linear 
piecewise  functions,  so  called  “the  roofs”  (Marchuk  and  Agoshkov,  1981;  Zienkiewicz,  1977). 

The  approximate  solution  is  represented  by  sums  of  the  form  •  Probe  functions 

i,k 

Kj,  i=l . 4,  TZj,  Jt5,  itp  are  also  chosen  among  .  The  vector  of  coefficients  will  be 

denoted  by 

It  is  obvious,  that  is  approximated  by  the  sum  of  (p/(A:,y)  . 

We  constantly  apply  the  so  called  “lumping-method”  to  derive  the  system  of  ordinary 
differentialequationsforvectors  ir,v,w,P,r,5,p,^.  IhsteadofGrammaliixGforfunctions  Kj^  j^  we 
use  the  diagonal  matrix  M  with  the  elements 

f  h 

Q* 

1,  ifi  =  / and k=m 
0,  ifi  ^lork^m 

The  Gram  matrix  Gi  for  functions  Tt^-^  in  a  space  with  the  scalar  product  weighted  by  /,•  is 
substituted  for  the  diagonal  matrix  C  with  the  elements 

(0(i,kUl,m)=^(i,kW,m)h  ■ 


where  = 


The  gradient  approximation  is  defined  by  the  formulae; 


M 


(Vip)a=  J  , 


a 


h  1=1 
KUi) 


<y\P)i,k-  J(I  • 

n*  >n=I  ^ 

The  trilinear  form  B  is  approximated  according  to  the  formula: 

M  M 


Qh  ^  Qh  m=/  l=I  ^ 

and  the  two  rest  terms  are  approximated  in  a  similar  way.  The  matrix  approximation  of  B  will  be 
denoted  as  A.  Applying  the  lumping  method  for  the  turbulent  mixing/diffusion  operators,  we 
obtain  the  matrices  D,  Dt  and  Ds  derived  with  the  use  of  the  approximate  equality 
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r 

dx  dx 

o!' 


M 


t=I 


dx  dx 


(and  the  same  for  the  last  terms).  The  coefficients  of  turbulent  mixing/diffusion  are  assumed  to 
be  constant  over  a  finite  element. 

As  a  result  we  have  the  following  system  of  ordinary  differential  equations: 

M  — +  Air-Cv=-— vJF  +  Dm  +  //  , 
dt  po 

f 

M  — +  Av  +  C«  =-— V*F  +  Dv+/2  , 
dt  Po  ^ 


v;p=sMp , 

-(V*)’'ff-(V5)’'v-(V5)''viT  =  /j  ,  (3) 

M^+Af=D,-f+/7-, 

dt 

M^+AS=D5S+/s  , 
p  =  p{T,S)  . 


1.1.3.  The  new  time  scheme 

In  the  previously  used  version  of  the  model,  there  was  a  time  step  limit  stipulated  by  the 
Coriolis  term  approximation  in  time,  and  there  was  the  explicit  “leap-frog”  scheme  for  the 
Coriolis  force.  According  to  the  CFL  criterion,  one  can  conclude  that  the  time  step  should  not  be 
greater  than  1.9  hours.  This  limitation  may  turn  out  to  be  too  severe,  therefore  the  next  step 
necessary  to  enhance  the  efficiency  of  the  model  is  to  change  the  Coriolis  term  time 
approximation.  We  used  the  Crank-Nicholson  time  approximation,  so  that 

CiT  = 

2 

The  new  scheme  leads  to  some  additional  difficulties  which  arise  because  the  matrix  for  sea 
level  was  nonsymmetric,  and  the  conjugate  gradients  method,  used  before,  is  no  longer  valid.  The 
nonsymmetricity  is  due  to  the  method  of  deriving  the  matrix.  To  obtain  the  system  of  equations 

for  ,  we  substitute  in  the  discrete  integral  continuity  equation: 

+  f  ^)dQ^  =  Fi ,  i  =  J . M 

\  dx  dy 
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which  is  derived  from  the  discrete  continuity  equation  by  means  of  the  summation  over  index  k  , 

KL(i) 

taking  into  account  the  equality  Vi  =  a  result,  we  have  a  system  of  linear 

k=l 


algebraic  equations  with  the  positively  defined  matrix  L  which  satisfies  the  equation: 

(Lo„, = J  <p^ao + g  J  iZ^i!k  I  + 


£2* 


1 

+- 

2 


dy 


, 

Bx 


o!'  ^ 


dy 


m=  I,..., 


M  . 


Here  we  use  the  denotations  a  =  2Ar  for  the  “leap-frog”  and  a  =  A/  for  the  Matsuno  time  step, 

^i,k  =  J  ’  h,k  =h^i,k->  ^i,k  =°i,k  +“  bij^  • 

£2* 

To  solve  the  problem,  we  tested  the  Sequence  Overrelaxation  method  (SOR),  commonly  used, 
and  rather  robust,  and  found,  that  for  a  time  step  greater  than  1.0  hour,  this  method  led  to  an 
underrelaxation  (it  converged  only  for  a  relaxation  parameter  of  less  than  1.0)  and  had  too  slow  a 
convergence.  If  the  time  step  is  larger  than  3.0  h,  the  SOR  does  not  converge  at  all.  For  a  time 
step  of  1 .0  h,  the  convergence  of  the  method  is  quite  acceptable  for  a  relaxation  parameter  of  1 .0. 


We  plan  to  examine  the  possibility  of  using  the  Generalized  Minimal  Residuals  method 
(GMRES)  with  preconditioning.  Where  the  matrix  for  sea  level  is  time  constant,  all  the  necessary 
calculations  of  the  preconditioning  matrix  and  Krylov’s  subspace  vectors  need  to  be  executed 
only  once.  We  hope  that  the  proper  choice  of  the  iteration  method  will  allow  us  to  use  time  steps 
of  about  6  hours. 


Note  2.  For  the  implicit  time  approximation  of  vertical  momentum  turbulent  mixing,  the 
matrix  for  sea  level  is  similar  -  only  the  right-hand  side  changes  slightly.  The  implicit 
approximation  may  be  very  useful  in  the  case  of  parameterization  by  Pacanowski  and  Philander 
(1981),  when  large  vertical  coefficients  (up  to  1&  m^s'  )  occur.  In  the  case  of  an  explicit  time 
scheme,  this  will  lead  to  extremely  small  time  steps  of  about  several  minutes. 


1.1.4.  Results  and  Problems 

As  a  test  of  the  new  numerical  scheme,  we  carried  out  two  sets  of  calculations.  The  first 
one  was  devoted  to  a  comparison  of  the  results  of  two  models  -  with  nonslip  boundary  conditions 
and  with  the  new  ones.  We  assumed  that  there  was  no  wind,  no  rivers  and  no  exchange  through 
the  passages.  Our  attention  was  focused  on  Fram  Strait,  rather  than  on  the  characteristics  of  the 
numerical  model.  The  results  are  shown  in  Fig.  1.1  and  Fig.  1.2.  They  show  that  our  efforts  were 
not  in  vain,  and  now  the  Atlantic  water  inflow  looks  more  realistic,  and,  in  addition,  it  can  be  seen 
that  the  overall  circulation  has  become  more  intensive. 
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Fig.  1.2.  Flow  velocity  at  250m  for  the  slip  condition. 


The  results  of  simulation  of  winter  circulation  specifying  the  wind  stress  and  the  water 
transport  through  the  passages  are  shown  in  Fig.  1.3  and  1.4.  Here  we  can  see  that  there  are  very 
“noisy"’  velocity  fields  at  the  boundary,  in  the  vicinity  of  relatively  intensive  currents,  especially 
off  Greenland  and  Scandinavia. 

Some  of  the  velocity  vectors  are  directed  opposite  to  the  main  currents  (in  the  region  of  the 
East  Greenland  current,  for  example).  This  effect  is  very  prominent  in  the  upper  layers  (compare 
with  Fig.  1.4  where  the  velocity  field  is  rather  smooth)  and  is  located  in  the  regions  of  strong  Joint 
Effect  of  Baroclinicity  and  Bottom  Topography  (JEBAR,  Sarkisyan,  1977).  Note  that  this  effect 
here  is  of  numerical  origin,  and  arises  from  the  inconsistency  of  the  bottom  topography  with  the 
temperature/salinity  fields  used  in  the  model.  The  T,  S  fields  are  rather  smooth  and,  in  the  region 
of  the  East  Greenland  current,  the  isotherms  and  the  isohalines  follow  almost  parallel  to  the  coast 
line.  The  bottom  topography  from  the  model,  here,  is  in  horizontal  steps  with  sharp  variations  of 
depth  from  point  to  point  and,  hence,  the  opposite  term 

/  " 

—  \j{H,p)dz 

Po  J 

(J  is  a  Jacobi  operator)  responsible  for  the  JEBAR,  is  unrealistically  large. 

We  also  examined  the  problem  of  velocity  oscillations  due  to  the  numerical  scheme.  Our 
finite  element  numerical  procedure  may  be  treated  for  the  case  of  uniform  triangulation  as  some 
finite-difference  scheme  on  the  A-grid  by  Arakawa’s  classification.  It  is  known  that  the  A-grid 
schemes  produce  numerical  noise.  We  may  expect  that,  where  the  boundary  condition  for  the 
velocity 

on  Fg  , 

is  approximated  in  the  horizontal  plane  by  three-point  correlation,  this  noise  will  be  strongest  in 
the  areas  of  sharp  coastline  variations  with  intensive  alongshore  currents.  To  examine  this 
assumption,  we  carried  out  numerical  experiments  with  the  minimal  sea  depth  being  limited  to  50 
m,  which  should  reduce  the  JEBAR  influence.  All  the  passages  were  blocked,  but  there  was 
baroclinicity  and  wind  stress  in  the  model.  The  result  is  presented  in  Fig.  1.5. 

Here  it  is  obvious  that  there  are  two-point  oscillations  at  the  side  boundary,  and  these 
oscillations  are  most  prominent  at  the  comer  points  of  the  model  area.  Finite-difference  schemes 
on  the  B  and  C  grids  do  not  produce  such  noise,  because  they  perform,  in  essence,  some  spatial 
averaging  over  boxes  (or  their  facets)  which  may  also  be  treated  as  cubic  finite  elements. 
Consequently,  we  may  consider  the  problem  of  a  noisy  velocity  field  to  be  one  of  interpretation  of 
the  results.  The  simplest  way  to  filter  noise  is  to  average  the  results  over  four  points.  The  result 
of  filtering  is  presented  in  Fig.  1.6. 
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Fig.  1.3.  Flow  velocity  at  the  25m  depth 
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Fig.  1.4.  Flow  velocity  at  the  500m  depth 
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Fig.  1.5.  Test  for  the  velocity  oscillations  at  the  side  boundary. 
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Fig.  1.6.  Flow  velocity  averaged  over  four  grid  points. 
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1.2.  Comments  on  the  choice  of  turbulent  mixing/diffusion  coefficients 

We  used  constant  coefficients  of  turbulent  mixing  and  diffusion  in  the  general  circulation 


model.  Their  values  were  as  follows: 

Horizontal  mixing  2T0  m‘s' 

Horizontal  diffusion  of  T,S  ITO  ’  m‘s‘' 

Vertical  mixing  50  m“s' 

Vertical  diffusion  of  T,S  1  m^s  ' 


One  may  argue  that  these  coefficients  are  very  large.  We  gathered  all  available  information 
on  the  coefficients  used  in  other  circulation  models  of  the  Arctic  Ocean.  The  results  of  the 
investigation  are  summarized  in  Table  1.1.  Because  of  poor  information  on  the  possible  values  of 
vertical  coefficients,  we  outline  only  the  horizontal  coefficients. 

Table  1.1.  The  coefficients  of  horizontal  mixing  in  different  models 


Model  and  horizontal 
resolutions 

Horizontal  mixing 

Recalculated  values 

Semjonov  G.,  Russia 

28  km  -  312  km 

Smagorinsky, 

8  10^-9  10 ’mV 

10®-  10^  mV 

Polyakov  I.,  Russia 

55.6  km 

Constant 

10  ’  m’s  ' 

410 ’mV 

Hakkinen  S.,  USA 

28  km  -  150  km 

Constant 

3T0’mV 

1.5  10 ’mV  -3  10® mV 

The  first  column  shows  the  designer  of  the  model  and  the  horizontal  spatial  resolution  used  in 
the  model.  The  second  column  shows  the  value  (or  the  range  of  values)  of  the  horizontal  mixing 
coefficient.  The  third  column  shows  the  effective  mixing  coefficient  recalculated  for  a  horizontal 
resolution  of  100  km.  It  is  not  very  easy  to  compare  these  data,  because  the  models  are  quite 
different  with  a  lot  of  specific  numerical  features.  In  Semjonov’s  model,  Smagorinsky’s 
parameterization  is  used.  Background  mixing  in  the  model  is  rather  small,  but  depends  heavily 
on  the  grid  size  and  the  velocity  deformation.  This  mixing  is  very  large  near  Greenland. 
Hakkinen’s  model  uses  very  smooth  fields  of  temperature  and  salinity,  and  even  in  this  case,  the 
mixing  is  close  to  the  results  of  our  model.  Polyakov  proposed,  formally,  the  less  viscous  model. 
However,  this  model  has  “numerical”  viscosity  due  to  the  first  order  schemes  applied  for 
advection.  Thus,  we  may  conclude  that  our  model  is  comparable  to  other  models  for  the  Arctic 
Ocean. 
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II.  CURRENT  VERSION  OF  THE  ARCTIC  BASIN  OCEANOGRAPHIC  DATA  BASE 


2.1.  Introduction 

This  section  describes  briefly  the  oceanographic  data  base  collected  to  support  the  thinking  on 
the  evolution  of  large-scale  anomalies  in  the  Arctic  Basin.  This  data  base  will  be  useful  for  many 
other  investigations  into  Arctic  Basin  oceanography. 

To  begin  with,  the  procedure  of  replacement  of  the  CTD  data  in  the  Arctic  part  of  the  World 
Ocean  Atlas  -  1994  (WOA-94)  is  described.  Then,  the  preliminary  results  of  combining  the 
Arctic  Master  Oceanographic  Observation  Data  Set  (MOODS)  with  the  WOA-94,  are  given,  and 
lastly,  miscellaneous  information  concerning  Arctic  oceanographic  data,  presently  available,  is 
listed. 


2.2.  The  new  version  of  CTD  data  in  the  Arctic  part  of  the  WOA-94. 

It  is  known  that  almost  one  third  of  all  the  WOA-94  CTD  profiles  were  truncated  due  to  a 
programming  error.  According  to  the  NODC  Internet  information: 

•  this  error  resulted  in  reporting  some  of  the  CTD  files  to  depths  shallower  than  they  were 
actually  measured; 

•  this  concerns  the  profiles  in  the  CTD  probe  files  only; 

•  the  error  has  been  corrected  and  new  CTD  files  are  currently  available  electronically  and 
on  CD-ROM; 

•  the  new  files  include  all  the  CTD  profiles  which  previously  contained  errors  as  well  as  all 
the  other  profiles; 

•  the  profile  identification  numbers  (PINs)  in  these  new  files  differ  from  the  PINs  on  the 
WOA-94  CD-ROMs. 

The  WOA-94  data  are  organized  by  ten-degree  latitude-longitude  squares  which  are  identical 
to  those  specified  in  the  ten-degree  square  numbering  scheme  of  the  World  Meteorological 
Organization  (WMO).  Within  each  of  the  WMO  subdirectories  there  is  a  file  containing  the  CTD 
profile  data.  Therefore,  the  ten-degree  squares  containing  CTD  files,  were  selected  from  the  CD- 
ROM  "Replacement  CTD  data"  (version  of  1996).  These  squares  were:  1800,  1801,  7800,  7801, 
7806,  7717,  7716,  7715,  7714,  7713,  1717.  Then,  all  the  CTD  data  located  in  the  corresponding 
squares  of  the  WOA-94,  were  chosen.  It  is  necessary  to  note,  that  not  all  of  the  squares  within  the 
Arctic  Basin  are  represented  in  the  new  version  of  the  CTD  files.  We  did  not  consider  the  squares 
corresponding  to  the  Nordic  seas,  the  Arctic  seas,  and  the  Canadian  Arctic  Archipelago  (7701- 
7710  and  1701-1715). 

To  characterize  the  differences  between  the  new  and  the  old  versions  of  the  CTD  profiles, 
three  types  of  profiles  were  grouped  (Table  2.1).  The  first  type  are  profiles  from  different 
versions  with  similar  country  codes,  expedition  codes,  dates,  times,  coordinates,  equal  or  more 
numbers  of  measured  characteristics  in  the  new  version,  and  with  more  numbers  of  measurement 
levels  and  deeper  levels  of  the  lowest  levels  of  observations,  in  the  new  version.  The  second  type 


IS 


are  profiles  from  different  versions  with  similar  country  codes,  expedition  codes,  dates,  times, 
coordinates,  lowest  depths,  and  number  of  measured  characteristics  and  levels,  while  the  third 
type  are  of  profiles  from  those  different  versions  which  are  not  included  in  groups  1  and  2. 

Table  2.1.  The  distribution  of  the  profiles  of  the  CD-ROM 
“Replacement  CTD  data  -  96”  and  the  WOA-94. 


Source  and  region 

Number 

of 

profiles 

Number 
of  profiles 
in  the 
first 
group 

Number 
of  PINs  in 
the  first 
group 

Number 
of  profiles 
in  the 
second 
group 

Number 
of  PINs  in 
the 

second 

group 

Number 
of  profiles 
in  the 
third 
group 

Number 
of  PINs  in 
the  third 
group 

1 

2 

3 

4 

5 

6 

7 

8 

Replacement  CTD 
data:  squares  1 800, 
1801.7800,  7801, 
7806 

155 

67 

67 

84 

84 

4 

4 

WOA-94;  squares 
1800,  1801,7800, 
7801,7806 

154 

67 

67 

84 

83 

4 

4 

Replacement  CTD 
data;  squares  -  1717, 
7713, 7714,  7715, 
7716.  7717. 

1653 

52 

51 

1545 

1542 

60 

60 

WOA-94;  squares  - 
1717, 7713,  7714, 
7715, 7716,  7717; 
*.CTD  files 

1653 

52 

52 

1545 

1541 

60 

60 

Analysis  of  Table  2.1,  and  further  examination,  allowed  us  to  conclude  that: 

•  the  new  version  of  CTD  profiles  replaces  the  files  with  name  *.CTD  in  the  WOA-94, 
rather  than  those  with  name  *.CTD2; 

•  in  addition  to  the  NODC  information,  we  have  found  other  differences  between  the  new 
and  the  old  profiles,  besides  the  number  of  levels  and  the  lowest  depth  of  observations 
(see  Table  2.1,  column  8); 

•  in  group  3  (column  8  of  Table  2.1),  the  observation  times  for  the  old  and  new  versions  are 
different  in  the  pairs  of  profiles  with  the  same  header  and  data; 

•  the  new  and  old  versions  of  the  CTD  files  contain  duplicates,  when  compared  by  country 
and  expedition  codes,  times,  dates,  coordinates  and  numbers  of  measured  characteristics 
(see  differences  between  columns  3  and  4,  5  and  6  in  Table  2.1),  while  there  are 
differences  in  the  data; 

•  about  20-25  %  of  the  PINs  are  the  same  in  both  versions,  but  correspond  to  different  CTD 
stations; 
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•  the  range  of  PINs  in  the  new  version  of  the  CTD  profiles  is  similar  to  that  in  the  WOA-94; 
therefore,  some  confusion  could  arise  if  one  were  to  substitute  these  CTD  profiles  into  the 
WOA-94,  using  the  new  PINs  for  the  comparison  of  the  headers; 

•  the  "Replacement  CTD  data"  CD-ROM  contains  one  more  station  than  the  WOA-94  (see 
Table  2.1,  column  2);  however,  the  temperature  and  salinity  in  this  profile  (PIN  101944) 
are  not  realistic. 

Fortunately,  the  other  errors  in  the  WOA-94,  detected  by  NODC,  and  listed  on  the  Internet,  do 
not  concern  the  Arctic  regions. 


23.  The  preliminary  results  from  combining  the  Arctic  Master  Oceanographic 
Observation  Data  Set  and  the  World  Ocean  Atlas  -1994 

There  are  9333  profiles  of  the  Arctic  Master  Oceanographic  Observation  Data  Set  (MOODS) 
in  the  public  domain.  The  World  Ocean  Atlas  -1994  (WOA-94)  contains  10366  profiles  within 
the  Arctic  Basin  margins  (15°W-105°E  80°N;  105°E-110°E  79°N;  110°W-130°W  77.5°N; 
130°W-140°W  78°N;  140°W-160°W  79.5°N;  160“W-170°W  77°N;  170°W-170°E  75°N; 
150°E-170°E  7rN;  135°E-150°E  70°N;  130°E-135°E  71°N;  125°E-130°E  75°N;  115°E- 
125°E  78°N;  90°E-115°E  80.5°N;  15°E-90°E  82.5°N).  The  general  information  on  field 
activity  in  the  Arctic,  has  allowed  us  to  assume  that  there  are  a  lot  of  duplicates  in  these  two  data 
collections.  Therefore,  the  comparison  between  the  MOODS  and  the  WOA-94  data  was  made 
before  the  merging  of  the  data  bases.  The  new  version  of  the  CTD  files  in  the  WOA-94  was  used 
for  the  comparison. 

The  common  parameters  which  compose  the  profile  header  in  both  the  MOODS  and  the 
WOA-94  are:  date  (including  day,  month,  year),  time,  coordinates,  number  of  measured 
characteristics,  and  number  of  observation  levels.  The  duplicates  were  determined  as  profiles 
with  similar  common  header  parameters,  except  for  time.  All  the  other  profiles  from  the 
MOODS  were  specified  as  new  for  the  joint  data  base  (Table  2.2). 


Table  2.2.  The  distributions  of  the  profiles  from  the  MOODS  and  the  WOA-94. 


MOODS  areas 

Number  of 
profiles  from 
MOODS 

Number  of 
profiles  from 
WOA-94 

Number  of  new 
profiles  from 
MOODS 

1 

2 

3 

4 

area  1;  70°N  -  80°N,  130°W  -  169°W. 

6409 

5388 

3754 

area  2,  80°N  -  90°N,  40°E  -  138°  W. 

1844 

3925 

955 

area  3,  80°N  -  90°N,  138°W  -  140°E 

111 

68 

47 

area  4,  70°N  -  80°N,  169°W  -  175°E 

764 

376 

534 

area  5, 80°N  -  90°N,  40°E  -  140°  E 

205 

108 

135 
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We  are  planning  to  apply  a  more  restrictive  algorithm  in  the  comparison  of  temperature  and 
salinity  values,  to  detect  additional  profiles  in  the  MOODS  which  do  not  coincide  with  the 
profiles  in  the  WOA-94. 


2.4.  The  current  content  of  the  Arctic  Basin  Data  Base 

Our  Arctic  Basin  Data  Base  is  growing  due  not  only  to  the  combining  of  commonly 
recognized  collections,  and  the  results  of  recent  expeditions,  but  digitizing  of  the  historical 
Russian  data,  also  provides  our  data  base  with  additional  information.  For  example,  we  have 
digitized  the  unpublished  archive  of  P.P.  Shirshov  in  the  framework  of  this  contract  (see  rows  17, 
26-41  of  Table  2.3). 

Table  2.3.  Miscellaneous  information  about  the  current  content  of  the  Arctic  Basin  Data  Base 


No. 

Data  description 

Location, 

deg.,  min.  or  region 

Date 

month., 

year 

Number  of 
real  stations 
or  modeled 
profiles 

mm 

2 

3 

4 

5 

1 

World  Ocean  Atlas  (WOA)  -  94 

Arctic  Basin 

06.1905- 

10.1990 

more  than 
10000  sts. 

2 

The  public  domain  of  the  Arctic  Master 
Oceanographic  Observation  Data  Set 
(MOODS) 

Arctic  Basin 

07.  1905- 
03.1995 

9333  sts. 

3 

Ice  Camp  “Turpan”  (Russia-USA) 

83.17N-83.35N 
23.45E  -  27.23E 

04.1994 

6  sts. 

■ 

Ice  Camp  “Narwhal”  (USA-Canada) 
(additional  data  from  “Narwhal”  are 
available  on  the  INTERNET). 

Lincoln  Sea, 
83.33N  -  83.35N, 
62.59W- 63.14  W 

04.1994 

19  sts. 

5 

"Oden-91",  Swedish  RV 
(only  small  part  of  observations) 

83.33N  -  83.35  N 
27.38E  -  27.58  E 

08.1991 

2  sts. 

6 

SCICEX-95,  American  submarine 
SSXCTD  data 

(SAIC  part  of  the  observation). 

section  between  72N, 
150Wand  85N,51E. 

04.1995 

63  sts. 

■ 

Winter  and  Summer  profiles  of 
Climatic  Fields  in  the 

Arctic  Ocean  (1994) 

Arctic  Basin 

3339 

profiles 

8 

“POLEX”  -  gridded  data  from  Russian 
quasi-synoptic  aircraft  surveys  of  the 
Arctic  Basin 

tracks  from  point 
Barrow  to  l)Svalbard, 

2) Franz  Josef  Land, 

3)  Severnaya  Zemlya 
and  from  the  Lincoln 
Sea  to  4)  Svalbard  and 

1  5)  Franz  Josef  Land 

03-05  of 
1955-56, 
1973-79 

711  profiles 

21 


9 


Ice  Camp  Gator  (the  exact  day  and 
time  of  measurements 
are  not  available) 

74N;  150W 

04.1992 

9  sts. 

Rose/Red  Ice  Camp  (the  exact  day  and 
time  of  measurements  is  not  available) 

73N;  164W 

04.  1986 

5  sts. 

Ice  Camp  SIMI  (no  salinity) 

74N;  150W 

04.  1994 

7  sts. 

Ice  Camp  FRAMIV  (only  SVPs) 

83N;  16E 

04. 1982 

5  sts. 

Ice  Camp  Tanya  (Approximate 
location,  time  and  date) 

86.5N,  14W 

04.  1986 

1  sts. 

Ice  Camp  Jerry  and  Shirley 

Lincoln  Sea 

04.1988 

16  sts. 

“Polarstem”  data 
(part  of  observations) 

81N-83N;30E-43E 

09.1983 

22  sts. 

“Polarstem”  data,  ARCTIS  IV/3 
expedition(part  of  observations) 

81.2N-  86.0N; 

20E  -  30E 

07  -08  1987 

18  sts. 

“North  Pole  -1”, 

Russian  drifting  station 

80.23N  -  88.53  N, 
06.15W-23.48E 

06.1937- 

12.1937 

28  sts. 

“North  Pole  -5”,  Russian  drifting 
station  (part  of  observations) 

85.24N-86.46N, 

82.23E-120.44E 

11.1955- 

11.1956. 

5  sts. 

“North  Pole-6”,  Russian  drifting 
station  (part  of  observations) 

82.25N-86.31N, 

53.40E-142.15E 

06.1958- 

03.1959. 

4  sts. 

“North  Pole-7”,  Russian  drifting 
station  (part  of  observations) 

85.41N-86.22N, 
44.26W-  157.25W 

01.1958- 

03.1959 

3  sts. 

“North  Pole- 13”,  Russian  drifting 
station  (main  part  of  observa-Ions) 

78.14N-85.55N, 
91.30E-  170.23E 

05.1965- 

04.1966. 

22  sts. 

“H  -  169”,  Russian  aircraft 
(part  of  observations) 

78.28N-81.37N 
176.44-  179.12  E 

04.1941 

2  sts. 

Some  results  of  Russian 
aircraft  expeditions 

80.04N  -  86.56N, 
157.20W-  142.09E 

04.-0.5  of 
1948,  1949, 
1950,  1955. 

16  sts. 

Russian  expedition  on  drifting  ice 
(small  part  of  observations). 

80.51  N  109.43  E 

04.1951 

1  sts. 

"Fram",  Norwegian  RV 
(famous  Nansen  expedition; 
part  of  observations) 

80.25N  -  85.30  N 
59.06E  -  132.04  E 

04.1894- 

11.1895 

4  sts. 

"Nautilus",  American  submarine 
(some  of  these  data  are  in  WOA-94) 

80.31N- 81.59  N, 
1.10E-24.45E 

08.  -  09 
1931 

9  sts. 

"Sedov",  Russian  RV 
(part  of  observations) 

75N  -  86N; 

4E  -  153.E 

08.1929- 

12.1939 

199  sts. 

"Litke",  Russian  RV 
(small  part  of  observations) 

81.40N,  9.23E 

81.1  IN,  94.55E 

09.1955 

09.1940 

mm 

"Stalin",  Russian  RV 
(part  of  observations) 

80.13N-81.00N; 

2E-4E 

12.1939- 

01.1940 

2  sts. 

30 

"Ermak",  Russian  RV 

79.07N  -  81.27N; 
4.40E  -  20.05E 

06-08  1897 

15  sts. 

31 

“Krasin”,  Russian  RV 

79.24N-  81.47  N; 
9.40E  -  45.00E 

07  -  09  1928 

31  sts. 

32 

“Maligin”,  Russian  RV 

81.37N-82.24N; 

51.54E-62.20E 

08  -  09  1932 

7  sts. 

33 

“Rusanov”,  Russian  RV 

79.15N-82.03N; 
54.00E-  101.32E 

08.1932- 

10.1936 

11  sts. 

34 

‘Taymir”,  Russian  RV 

08.-09.  1932 

15  sts. 

35 

“Persey”,  Russian  RV 

80.00N  -  80.0 IN; 
3.03E-  10.35  E 

08.1934 

4  sts. 

36 

“Sadko”,  Russian  RV 

77.00N-  82.41N;3.19E 
-  94.38E 

08.1934- 

09.1936 

101  sts. 

37 

“Sibirakov”,  Russian  RV 

79.16N- 80.17  N; 
6.40E-91.34E 

09.1936- 

06.1939 

5  sts. 

38 

“Nerpa”,  Russian  RV 

79.00N-82.15N; 
50.23E  -  95.23E 

08.1937- 

09.1940 

44  sts. 

39 

“Belgica”,  RV 

79.12N-80.17N; 
1.29E-  17.02E 

06.  -  07.  1905 

13  sts. 

40 

“Farm”,  RV 

79.10N- 80.16  N; 
8.10E-  12.05E 

07.-08.  1910 

10  sts. 

41 

“Veslemov”,  ?  RV 

79.05N  -  80.39N; 
V.OIE-  17.14E 

07.-08.  1912 

39  sts. 

The  Data  Base  described  above,  will  be  used  for  analysis  of  the  evolution  of  large-scale 
anomalies  in  the  Arctic  Basin. 
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